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Abstract. We analyze two configurations for laser cooling of neutral atoms whose 
internal states store qubits. The atoms are trapped in an optical lattice which is placed 
inside a cavity. We show that the coupling of the atoms to the damped cavity mode 
can provide a mechanism which leads to cooling of the motion without destroying the 
quantum information. 
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1. Introduction 

Quantum optical scliemes for quantum information processing employ long lived internal 
(spin or hyperfine) states of atoms or ions to store the qubit [J- In many of the proposals 
for scalable quantum computing these qubits are stored in a storage area, and must be 
moved together to perform e.g. a two-qubit gate operation entangling the atoms |21IS1- 
Moving atoms will typically cause heating, and thus the question arises to what extent 
the motion of the qubit atoms can be cooled nondestructively, i.e. without affecting the 
coherence of the qubit. In the case of quantum computing with trapped ions this is 
achieved by sympathetic cooling the "qubit ion" is coupled via the Coulomb force 
to a second "cooling ion" (of a different species) , which is laser cooled [21 E] • This cooling 
is nondestructive because the Coulomb force is identical for the two qubit states |0), 
Quantum computing with neutral atoms can be performed by storing atoms in optical 
lattices generated by counterpropagating laser beams. By moving the atoms in a spin 
dependent optical lattice [3 |H] one can induce controlled collisions between the qubits 
leading to multi particle entanglement [71 HO]- One possibility for nondestructive 
cooling in an optical lattice is sympathetic collisional cooling by immersion of the qubit 
atoms stored in the lattice in an atomic Bose Einstein condensate ^I]. Of course, laser 
cooling of qubit atoms stored in the wells of an optical lattice potential cannot be directly 
applied, as the optical pumping associated with the laser cooling process destroys the 
coherence. Instead we will investigate below a scenario of nondestructive laser cooling 
of atomic qubits in moving optical lattice potentials which is achieved by coupling to a 
(low-Q) optical cavity. These results are also of interest for implementations of quantum 
computation and quantum communication with optical cavities |12[ IT^ [Til El dj • 

Cavity assisted laser cooling is a well-studied process [T7 1 ITH l ITH l 1 ^ 1 ^ 1 ^ 1^ . 
both theoretically and experimentally. In the present paper we will elaborate on two 
distinguishing features which make cavity assisted laser cooling of interest for quantum 
computing with atoms in optical lattices. We will study a scenario where a set of qubit 
atoms is stored in individual wells of an optical lattice potential. By moving the lattice 
we will heat the atom motion. However, to the extent the optical lattice potentials seen 
by the two qubit states are identical, this transport will not destroy the coherence of 
the internal atomic states. Cavity assisted recooling of the motional degrees works by 
converting phonons to cavity photons in a coherent laser assisted process. These cavity 
photons dissipate through the cavity mirror. With an appropriate choice of internal 
qubit states, this cooling process does not distinguish between the two qubits, and is 
thus nondestructive. We also note that moving N qubits in an accelerating lattice will 
excite the collective center of mass mode of these atoms. This is precisely the mode 
which is cooled most efficiently by our cavity cooling scheme. The scheme we investigate 
is, in principle, capable to achieve ground state cooling in the optical wells, although in 
practice nondestructive cooling to the Lamb-Dicke regime, where the atom is localized 
in a region much smaller than the optical wave length, may be sufficient. 

The paper is organized as follows. In Sec. El we present two different specific setups 
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and we calculate the equations of motion in both cases for a single qubit and for the 
center of mass motion of N qubits. In Sec. El we derive analytical expressions for the 
steady state temperature Tf and the cooling time r of the atoms and compare them 
with numerical calculations. Finally, in Sec. |3] we discuss possible imperfections of the 
system. 

2. The Model 

In this section we introduce two different setups for cooling the motion of neutral atoms 
without destroying the quantum information stored in two long lived degenerate internal 
ground states \gi) = |0) and \g2) = In the first setup the atoms are trapped and 
cooled by the field of a ring cavity while in the second case trapping is achieved by an 
optical lattice and the cavity is solely used for cooling the atoms. To be specific, we 
will choose the two qubit states to be two Zeeman states, e.g. \ J = 1/2, M = ±1/2), 
and we will assume identical (symmetric) coupling of these states to the laser and the 
cavity mode to achieve cooling which does not distinguish between the two qubit states. 
Our goal is to derive cooling equations first for a single qubit, and then for the center 
of mass motion of qubits. 

2.1. Optical Lattice in a Ring Cavity 

We consider neutral atoms in a ring cavity as shown in Fig. ^). Two coherent fields 
Pr and Pl drive two counter-propagating running waves adding up to a standing wave 
inside the cavity. We first consider a single qubit in this ring cavity and then extend 
the system to N qubits. 

2.1.1. Single Qubit In a frame rotating with the laser frequency ujl, the Hamiltonian 
describing the cavity fields and one atom is given by 

H = Ha + Hc + Hac + H-p + Hqf + -f^bath + -f^sys-bath (1) 

with {h = 1) 

HA = ^ + {-A){\e,){e,\ + \e,){e2\), (2) 
He = (-AJ(aJja/j + a^ai), 

Hac = 9 [{ane'^"''--^'^ + aie-'('=^+-^*)) (a+ + a+) + h.c] , 

j=R,L 

HcF = iv^ 5^ [diu [a%,,e'^^' - h.c.) . 

j=R,L ^ 

Here the detuning of the laser from the atomic transition frequency uj^g is given by 
A = ul — and the detuning from the cavity frequency Uc by Ac = uji — Uc- 
Furthermore, Ha gives the kinetic and the internal energy of the atom, respectively. 
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where p is the momentum operator and m is the mass of the atom. The two long 
hved degenerate ground states carrying the qubit are \gi) = |0) and \g2) = |1), as 
sketched in Fig. [T] and we assume the atom to have two degenerate excited states |ei) 
and 1 62). The Hamiltonian Hq describes a ring cavity with a clockwise and counter- 
clockwise propagating mode, and bosonic destruction (creation) operators aR(^L) ('^r(l))- 
In dipole and rotating wave approximation, the interaction of the atom with the cavity 
field has the form Hp^c, where = \ei) {gi\, crj = |e2) ((72! and aj = {ctj^Y , j = 1,2. 
Furthermore, g denotes the single photon Rabi frequency, and k is the wave vector of 
the laser light. The free Hamiltonian for the external driving field is denoted by Hp 
and Hqf describes the interaction of this field with the two cavity modes, where k is 
the photon loss rate through the cavity mirrors. The bosonic annihilation (creation) 
operators for the field in the mode u, b^^Fi.{L) i^lj r(l)) fulfill the commutation relations 



[buj,R{L), = ^(^ ~ ^')- Finally, the term -ffgys-bath describes the interaction of the 

atom with a photonic bath ifbath, which leads to spontaneous emission. These terms 
are of standard form and can be found e.g. in I25j. 

We assume the electric field driving the two cavity modes to be coherent, 
representing a classical field with amplitudes An,R(L)(^)- As described in detail in 
[Appendix A we go to a picture where the incident driving field explicitly enters the 
Hamiltonian. This results in a vacuum input state for the cavity and the additional 
driving term 



CF 



2k {l3,^{t)Al - A+An(t)* 



(3) 



in the Hamiltonian. Here we have defined A± = (a^ ± a^) / V2, and furthermore 
assumed symmetric driving Pin,R(t) = Pm,L(t) = Pin{t). 



a) 



b) 



/3l 




A 



A, 



Figure 1. a) Experimental setup as described in the text, b) Schematic picture of 
the internal level structure, the involved frequencies and detunings. The ground state 
Iffi) (I52)) is identified with the qubit state |0) (|1)). As a concrete example one can 
chose the fine structure states (e.g. of Na or Rb) as shown in Fig. 



The dynamics of the system thus obeys a master equation (see e.g. 
p = -i [Hsys,p] +C{p), 



(4) 
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where the Hamiltonian of the system is given by H^ys = Ha + Hq + i^Ac + -f^cF- The 
Liouville superoperator C has the form 

^(/^) = f E - Map - p4 a) (5) 

i=± 

+ I E - ^t^iP - P^>i) ■ (6) 

i=± 

The first hne Eq. ^ describes loss due to the imperfectness of the cavity mirrors, and 
the second hne Eq. is due to the spontaneous emission of the atom with emission 
rate 7. Note that based on the choice of atomic states and laser and cavity polarizations, 
the laser and cavity couphngs as well as the spontaneous emission rates are assumed to 
be identical for the two qubit states. 

Next we decompose the cavity mode operators A± into a C-number a± = {A±) 
describing the coherent part and an operator part 6A± = A — a describing the 
fluctuations. The mean values {A±) = a± are calculated in the steady state. From 
Eq. (jH) the respective equations of motion for {A±) can be easily obtained, yielding 

= (iA, - ^) + V^An, 

A(A_) = (iAe-^)(A_). (7) 
From the steady state solution of Eqs. ((Zj) we find 

a+ = — = a, a_ = 0. (8) 

-lA, + 1 ^ ^ 

Applying this transformation to the system Hamiltonian we find ifgys = Hj^ + Hc + H[^t, 
where in He the operators A± are replaced by SA± and i^int is given by 

Hint = V2g [{{6A+ + a) cos kx + iM_ sin kx) (ctj*" + cr^) 
+h.c.] . 

The Liouville operator C reduces to 

^(/^) = f E (25^iPMl - SAUAp - p6Al6A,) , (9) 

where we have assumed large detunings A and consequently neglected the spontanous 
emission, as it is suppressed by a factor oc g'^/A'^. Furthermore, in this approximation 
we can adiabatically eliminate the excited states |ei) and |e2) and obtain 

/ 2 

Hes=h®(^-A^ (^5Al5A+ + mLM_) (10) 



2g^ 

+ 



(^a^ cos^ kx + ia{6Al_ — SA_) cos kx sin kx + 5A^5A_ sin^ kx^ ^ 



A 

where we have only kept the leading terms in 6A, assuming small fluctuations | (6 A) \ -C 
In Eq. (ITn|l the operator Ij = \gi) {gi\ + \g2) {g2\ which we do not write explicitly 
in the following. Thus, due to our choice of symmetric coupling and dissipation, the 
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information encoded in the qubit (i.e. the amphtudes of the two internal ground states) 
is not affected by the time evolution Eq. Furthermore, the operators 6A^ and 6A\_ 
do not couple to any other degrees of freedom of the system and, therefore, we drop 
them in the following. 

The periodic optical lattice inside the cavity is given by the term Vq cos^ kx in 
Eq. (fnH) with a depth of Vq = 2g'^a'^/\A\. For a blue (red) detuned lattice and 
a kinetic energy of the atom much smaller than the depth of the optical lattice 
(i.e. (p^) /2m ^ Vq), the atoms are trapped near the (anti)nodes of the standing wave. 
Thus, we can approximate each well of the lattice by a harmonic oscillator with frequency 
v = 2gak/\/mA. In this approximation C is given by 

C{p) = — (2apa^ — a}ap — pa}a) , (11) 

where we have introduced a = (i4_ and the system Hamiltonian reduces to 

-f^sys = + VcC^a + i(yfefr (c + c^) {a} — a) . (12) 

The bosonic creation (annihilation) operator (c) describes the quantized motion of 
the atom in a harmonic potential of frequency v coupled to the cavity mode a with the 
effective couphng strength g^-s = v /2ar}, as shown in Fig.|21 In Eq. (fT^ we also used the 
definitions Vc = —^c + ^ for the frequency of the cavity mode a, and rj = ■sJkP' /2mi' 
the Lamb-Dicke parameter. The shift 2g^ j A can be interpreted as a refractive index 
due to the back-action of the atom on the cavity field. In the limit (yf/z/ -C 1, i/ = z/c we 
can apply the rotating wave approximation (RWA) to Eq. (HH), and find 

Hbwa = i^c^c + VcO^^a + i^feff (ca^ — ac^) . (13) 




Figure 2. The two coupled harmonic oscillators corresponding to the atomic motion 
and the cavity with frequencies u and Uc, respectively. The coupling strength is denoted 
by ffcff and k, is the photon loss rate through the cavity mirrors. 

We finally turn to the physical interpretation of our results. The coupled oscillator 
Hamiltonian Eq. (fT^ describes the coherent transfer of phonons (vibrational excitations) 
to cavity photons. Due to its high frequency, the optical cavity mode is coupled to an 
effective zero temperature heat bath. Cavity phonons will decay through the mirror 
according to Eq. (fTTj) . dissipating the energy. 

A physically equivalent picture can be given in terms of a dissipative optical lattice 
according to Eq. (fTUI) . In the ring cavity we have two counterpropagating modes, which 
we rewrite as cos /ex (+) and sin /ex (-) standing wave modes. By driving the cavity 
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symmetrically from the outside with two lasers, the standing wave coskx will have a 
large occupation of phonons jap, giving rise to an optical lattice potential, while the not- 
driven sin kx mode will be in a vacuum state. Motion of an atom in the optical lattice 
will redistribute the photons among the two modes. In particular, atomic motion will 
transfer photons to the (initially empty) sin /ex mode, i.e. the optical lattice potential 
will be slightly shifted by atomic motion. On the other hand, the photons occupying 
the sin kx mode can decay through the cavity mirror, giving rise to a damping. Thus 
we have the picture, where the atomic motion in the optical wells acts back on the 
lattice to produce a (small amplitude) displacement of the lattice. These oscillations 
of the optical lattice will damp to their equilibrium position by the cavity decay, thus 
cooling the motion of the atom. With an appropriate choice of qubit states we assure 
that the dissipative optical lattice potential seen by the two qubit states are identical. 
This ensures that the cooling does not destroy the coherence of the qubit. 

Finally, we note that our laser cooling model maps on to a remarkably simple 
mathematical model of coupled harmonic oscillators with linear damping. Thus we find 
a straightforward solution in terms of Gaussian states of the two modes. Below we will 
extend the model to N qubits, and then to the case of a time dependent moving lattice 
to illustrate the cooling dynamics. 

2.1.2. N Qubits The derivation of the A^-particle Hamiltonian is similar to the 
calculations presented in Sec. 12.1.11 and leads to 

N N 

H = i^Yl + '^c(A^)a^a + i^efr ^ [c^ + cj.) {a} - a) , (14) 
/x=i M=i 

where the index /i labels the particles. All parameters apart from udN) = — + 
2g'^N/ A and the refractive index now given by 2g'^N/ A remain unchanged. In Eq. (fT^ 
we have assumed the effective coupling constant to be the same for each particle 
igeS,ti = des)- The Hamiltonian Eq. (fT^ describes N harmonic oscillators corresponding 
to the quantized motion of the atoms coupled to the single cavity mode a. The coupling 
term in Eq. (|Tl|l has exactly the same form as in the single qubit case Eq. (|T2|) if 
we define the center of mass (COM) annihilation (creation) operator C (C^) given by 
C = ^/VN ^fj^Cf^ and introduce the coupling constant g = y/Ng^s- The resulting form 
of Eq. (fT^ already indicates that only the collective mode C will be cooled in our 
approximative description. We note that previous pubhcations [211 12Z| have already 
shown in a semiclassical approximation, that also the (N-1) relative vibrational modes 
can be cooled with our setup, but on a much larger timescale. 

As will be discussed in detail below, transport of atoms in an accelerating lattice 
leads to coherent oscillations of the center-of-mass mode of all the atoms. Our cooling 
scheme cools this particular collective mode with high efficiency due to the collective 
enhancement ges — > gesV^ of the phonon-photon. 
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2.2. Optical Lattice in a Standing Wave Cavity 

As a possibly more convenient setup, we replace the ring cavity of Sec. l2.1l bv a standing 
wave cavity providing the cooling. We assume the optical lattice holding to the atoms to 
be produced by two counter-propagating laser beams, which we shine into the resonator 
at a small angle, as shown in Fig. El whereas in the previous system the driving field 
entered the cavity through the cavity mirrors. The total Hamiltonian takes on the form 




Figure 3. Schematic experimental setup as described in the text. 

H = Ha + Hc + Hac + ^al, where 

^A = |^ + (-A)(|ei)(ei| + |e2)(e2|), (15) 

He = (-AJa^a, 

Hal = ^i sin kx (cr^ — cr~) , 

Hac = 9 cos kcX (^aa~^ + aV~) , 

and 0"^ = {af + o"^). The Hamiltonians Ha and He denote the energy of the atom and 
the cavity with bosonic creation and annihilation operators and a, respectively. The 
interaction of the laser with Rabi frequency Q and the atom is denoted by Hal and i^Ac 
describes the atom-cavity coupling {g again denotes the single photon Rabi frequency). 

Applying the same approximations as in Sec. 12.11 we find an identical model with 
the replacements i/ — z/ = Q'^rj'^/A and Uc ^ I'c = —Ac + fi'^/A, as well as the 
effective coupling g^s (jcs = g^ri/2A in Eq. (jl)). The Lamb-Dicke parameter rj is 
rj = yjkj^mv ~ kcjlmv.^ where k^ and k denote the wave vectors of the cavity and 
the laser light, respectively. A straightforward extension to N qubits results again in a 
Hamiltonian of the form (|THl . 

The physical interpretation of this setup is very similar to the previous one with 
the only difference that the optical lattice is formed by two additional lasers instead 
of the cavity itself. In this case the motion of the atoms transfers photons from the 
laser field into the (initially empty) cavity. As in the previous setup the decay of the 
cavity photons through the cavity mirrors leads to cooling of the atomic motion without 
affecting the qubit state. 
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2.2.1. Accelerated optical lattice Finally, we generalize the model to an accelerated 
optical lattice. In the setup of Fig. El the lattice can easily be moved by introducing a 
time dependent relative phase (j){t) between the two counter-propagating laser beams 
[3 IH] . We find a system Hamiltonian 

H = vc^c + z/cfl^a ^ (c + c^) — iga cos {4>{t)) (a^ — a) 

+ ig cos (0(t)) (c + c^) {a^ - a) + —c?-, (16) 
with z/c = — Ac + (7^ cos^ {(p)/ A, a = (p/rj. The values for u, g and t] are defined as above. 



3. Results 

Below we solve the coohng equations Eq. (jH). In Sec. l3.1l we derive analytical expressions 
for the atomic steady state temperature Tf and the cooling time r. In Sec. IH.2l we show 
examples for the numerical solution of Eq. Q in the case of a static lattice as well as 
for an accelerated lattice. 



3.1. Analytical Calculations 

The master equation ^ consists of a quadratic system Hamiltonian and a linear 
damping term. Thus, the Heisenberg equations for the expectation values of the second 
order moments can be found easily (see Appendix B[ Eq. ()B.1|) ). From these we calculate 



the atomic steady state temperature Tp 

kBTf = Eo = u{c^c)^^= (17) 



g!^ , + 



2u ^ 16z/e ^ u [k^u + 4z/e(-4^efr + ^^c)] ' 
with kB the Boltzmann constant and Eq the steady state energy. We note that the 
system is always cooled to the ground state in the limit g^s -C z/ ~ z/g (i.e. where the 
RWA Eq. (O is vahd). 

The atomic energy -E'atom(^) = {^^^)t '^^^ written in the form 

^atom(t) =5^Cfce^^* + Eo, (18) 
k 

where the are eigenvalues (all having negative real parts) obtained from the equations 
of motion and the coefficients are found from the initial motional state of the atom. If 
one term dominates in the sum Eq. p8|) then l/lReAfc] with the corresponding eigenvalue 
Afc gives the cooling time r. In general, however, the cooling time will be determined 
by the eigenvalue A^ with the largest real part unless the corresponding coefficient 
vanishes. We therefore find an upper bound for the cooling time r by 

r = . , |. • (19) 
mm||ReAfc||fc 

In RWA we can easily calculate the eigenvalues A^ from the equations of motion 
given in Appendix B[ Eq. ()B.2|) and find simple analytic estimate for the cooling time 
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r in two interesting limits. In the "Doppler limit", i.e. ^ ^ k, z/^ = I we find 
T ^ k/AqIq and fc^Tj ^ -I, in close analogy to free space Doppler cooling (see e.g. |2H1)- 
In the "sideband limit", i.e. ~ k -C z/ ~ z/c, we obtain t ~ ^ and ksTf ^ §1^/21'. In 
the numerical examples in Sec. 13.21 one can see that those analytical expressions are in 
very good agreement with the exact numerical calculations. These results for the single 
qubit case can be directly applied to the COM motion of N qubits by replacing g ^ g. 

3.2. Numerical Results 

In this section we show the numerical solution of Eq. for different parameter regimes 
in the case of a static lattice as well as for an accelerated lattice. Fig. HI shows the 



a) 6 




25 50 



Figure 4. a) Energy of the atomic motion i?atom and b) energy of the cavity mode 
^'cavity as a function of time vt. The parameters are chosen K = l,i/ = i^c = 2 and 
5eff — 1/8. The dashed Unes show the time evohition with the effective Hamihonian 
i?eff in RWA, the sohd Hnes are without the RWA. 

time evolution of the atomic energy -E'atom(i) and the energy of the cavity mode a, 
-^'cavity = ^s a fuuctiou of time vt. Initially the system is assumed to be in a 

product state with the cavity mode a in its ground state and the mode for the atomic 
motion in a coherent state. The time evolution of the two energies shows that energy 
is transferred from the mode c for the atomic motion into the cavity mode. The energy 
of the cavity mode decreases due to the photon loss through the cavity mirrors, which 
leads to a cooling of the atomic motion as well. Furthermore Fig. |3] shows a comparison 
between the time evolution in RWA and without RWA. We can see, that in the RWA 
no oscillations of the energies occur, i.e. in this approximation the energy is transferred 
smoothly from one mode into the other and dissipates at a rate k. 

In Fig. we compare the analytical results for the cooling rates and the final 
temperatures with the exact numerical solution in two typical examples. One can see 
that the cooling times which we have estimated by using Eq. ()13j) . fit the exact numerical 
calculations very well, and the final state temperature (which is almost zero in Fig. Et- 
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Figure 5. Energy of the atomic motion i^atom/i^ versus time I't. The parameters are 
chosen a) k = O.l, v — Vf. = QcS = 0.1 and h) n — 1, v — ly^ = 1 and g^g = 0.1. The 
dashed curves show an exponential decay with using the analytical expression for r. 
The dotted line in b) gives the analytical value for the final energy Eq/v. 



coincides exactly with the analytical value calculated from Eq. (jl7j) . 

a) 10 f 
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Figure 6. a) Energy of the atomic motion Ei^tom/v versus time vt. Initially the 
atom is in its motional ground state, then motional excitations are created by non 
adiabatic acceleration of the lattice. Afterwards cooling is switched on at vt = 51 with 
parameters n ^ \, v = Vc = ^ and gcff = 0.1. b) Expectation value x{t) — {x{t)) of the 
atomic position in dimensionless units and c) the relative phase between the two 
lasers. 



Finally, Fig. IHl shows the excitation of the motional state of an atom, initially 
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assumed to be in its motional ground state, by non adiabatic acceleration of the lattice. 
The chosen velocity profile is typical for a two qubit quantum gate where one atom is 
moved close to a second one for performing the gate operation. While here we assume 
that the motional state is excited during one such movement it might typically only be 
necessary to re-cool the atom after a few gate operations have been performed. After 
the motion has finished we turn on the cooling to bring the atom back into its motional 
ground state. The results are summarized in Fig. IHl The time evolution of the atomic 
energy 

Ea = V (etc) - — ((c) + (ct)) + -a^ - — ((ct) - (c)) + — , (20) 

which includes the kinetic energy part due to the motion as well as the vibrational 
excitation inside the trap is shown in Fig. IHK- The motion of the optical lattice is 
caused by a time dependent relative phase 0(t) = kxo{t) (c.f. Fig. with xo(t) the 
mean equilibrium position of the atom. During the acceleration, the mean atomic 
position x{t) = {x{t)), shown in FiglHfc, follows its equilibrium position Xo{t) quite well. 
After the acceleration has stopped and the cooling has been turned on the motion of 
the atom is cooled down to the temperature Tj (EqlTTj) within a few trap periods. 

4. Imperfections 

In this section we discuss sources of imperfections leading to qubit errors for a specific 
internal level structure of the neutral atoms storing the qubits. We consider the fine 
structure levels of an alkali atom as e.g. Na or Rb shown in Fig. [7| and the laser light 
ideally travelling along the quantization axis z. We assume the qubit state |0) (|1)) to 
be encoded in the nij = —1/2 {mj = +1/2) state of the 5*1 manifold. This qubit state 
is coupled to the excited state mj = —3/2 {nij = +3/2) of the Pa manifold by a a~ 
((T+) polarized laser with detuning Ai (A2). We consider two kinds of imperfections 
(i) qubit flip errors due to misalignment of the laser leading to small errors in the 
polarization vector 6e and (ii) stray magnetic fields 5B parallel to the quantization axis 
or equivalently imbalances in the two Rabi frequencies leading to qubit phase errors, 
i.e. the accumulation of a relative phase between |0) and |1) during the cooling process. 

The probability Pfjip for a spin flip error (|0) ^ |1)) to occur in a given time interval 
r is found to be 

= T^-'^^3, (21) 

where Se^ is the component of the polarization vector in the direction of the quantization 
axis. Inserting typical values r ~ 20/z/ (cooling time found for the parameters used in 
Fig. HI) and 863 ~ 10"^ we find Pmp ~ 10"^ 

A magnetic field in z-direction SB3 leads to a Zeeman shift SE = gmjfiB^B, with 
fiB denoting the Bohr magneton and g the Lande factor of the different states. The 
resulting difference in the transition frequencies from |0) ^ |ei) and |1) <-> |e2) leads 
to two different Stark shifts f2^/Ai and Q'^/A2. Note that imbalances in the Rabi 
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frequencies or differences in the trapping potentials of the two qubit states also lead to 
differences in these Stark shifts and can thus be treated in an identical way. The relative 
phase accumulated in time r is given by 



^rel 



2A 



TfisSBs. (22) 



^2 

For the typical values used above and 6B-^ ~ IG we obtain a relative phase (pj-ei 
5 X 10-'\ 





2 

Figure 7. Fine structure and laser configuration as described in the text. 



5. Conclusion 



In this paper we have studied the nondestructive cooling of atomic qubits in an optical 
lattice by coupling to a dissipative optical cavity. In particular, we have proposed 
two specific experimental configurations, where the optical lattice was produced in the 
interior of an optical cavity, and the required dissipative mechanism was provided by 
photon loss through the mirrors of the cavity. In the first setup two counter-propagating 
laser beams were driving two modes of a ring cavity, resulting in an intra-cavity standing 
wave. We interpreted the resulting cooling in terms of a dissipative optical lattice which 
by an appropriate choice of atomic states became independent of the qubit state. In 
the second setup the ring cavity was replaced by a standing wave cavity and an optical 
potential was generated by two counter-propagating external laser beams intersecting 
the cavity at a small angle. In both cases the physics of cooling could be described 
as two coupled oscillators consisting of the phonon mode coupled parametrically to the 
optical cavity mode, which was damped by cavity loss providing the dissipative element 
for the cooling process. 

Transport of atoms in the lattice gives rise to coherent oscillations of the center- 
of-mass mode of all the atoms. Our cooling scheme cools this particular collective mode 
with high efficiency due to the collective enhancement g^s gcsVN of the coupling 
of the phonon to the photon modes. We find cooling times which are similar to those 
observed for other schemes ^T] and are much smaller than the relevant decoherence 
times for the internal atomic states used in quantum information processing with 
neutral atoms. Therefore our scheme can be applied repeatedly to cool the motion 
of atoms which was excited during gate operations that involved moving the atoms. 
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While the present work has emphasized apphcation of this coohng scheme in a quantum 
information context, couphng the motion of atoms in optical lattices to photon modes 
of a cavity might also be interesting from the perspective of quantum degenerate gases 
in optical lattices. 
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Appendix A. Transformation to a Classical Driving Field 

We describe the transformation which allows the incident field to be represented by 
a classical driving field in the Hamiltonian (^, similar to the work by Mollow |29j . 
The state vector \ip{t)) giving the system and the driving field satisfies the Schrodinger 
equation 

i^^im) = H\m), (A.i) 

where H is given by Eq. (^Q). The initital state is assumed to be a product state of the 
system and a coherent field state 

1^(0)) = l^sys(O)) DiiP^^n, Pu.,l}) |vac) , (A.2) 
where |vac) is the vacuum state and 

mPu^^R, = exp ( 5^ fdu (^P^.bl^ - ) (A.3) 

\j=R,L J 

is the coherent displacement operator (see e.g. ppj). In order to transform the coherent 
incident field to a classical driving field we apply the unitary transformation 

im) ^^({/?.,/?e--*, /?.,z.e-^"*}) |^(t)) , (A.4) 

which results in an initial vacuum state for the driving field and in the transformation 
K,R{L) K,R(L) + /5a;,/?(L)e~"^*. Assumiug Symmetric driving (3^^^r = (3^^l = Puj the 
Hamiltonian now has the form Eq. (P) with the additional term Hqf defined in Eq. (jH)) 
and we have used 

An(t) = / da;e^('^^--)*/3^. (A.5) 
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The Heisenberg equations for the expectation values of the second order moments can be 
obtained from the master equation Q with the system Hamiltonian Eq. p2|) . yielding 

^(c^c) = {{ca) - (ca^) - (c^a) + (c^a^)) , (B.l) 
^{a^a) = 9cfi {{ca) + (ca^) + (c^a) + (c'^a"^)) , 

— (ca) = gcs ((c^c) + (a^a) + (c^) - (a^)) + (-iz/ - iz/^ - -) (ca) + ges, 



(ca^) = g^s ((c^c) - (a^a) + (c^) + (^a^"'^) + (-iz/ + iz/, - |) (ca^) 



dt 
d_ 

dt 
d 

dt 
d_ 

dt 

^(c^) = 2g,s (- (ca) + (cat)) + 2iz/ (c^) 



— (c^a) = g^s y{c^c) - (a^a) + (^c^ ^ - (a^) j + (iz/ - iz/, - -) (c^a) , 

(c^at) = g^s (^(c^c) + (a^a) + ^c^'^ - ^a^'^) + (iz/ + iz/, - |) (c^a^) + c,, 



dt 
d / 



(a^) = 2c/eff ((ca) + (c^a)) + (-2iz/, - k) (a^) 



^(ct^)=2^efr((cta>-(cat»+2iz/(ct^), 

(^at'^ = 2g,s ((c^at) + (ca^)) + {2iu, - zt) (^a^'^ . 



dt 



In RWA the time evolution of the expectation values of the quadratic operators is 
given by a homogeneous set of linear differential equations, which can be written in the 
form 



-y{t)=Mm- 
Here y{t) is given by 

y(t) = ((c^c) , (a^a) , (cat) ^ (^c^a))^ . 
and the time evolution matrix Ai has the form 



M 



(B.2) 
(B.3) 



/ 





-5'cfr 


-5'efr 







— K 


5'efr 






Qcs 


-fi'eflf 


i(z/c - z/) - f 







\ 9cs 


-fi'eff 





i(z/ - z/J - 


f / 
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